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SECTION  1 
INTRODUCTION 

This  study  was  conducted  for  the  Air  Force  Materials  Lab¬ 
oratory  to  analytically  determine  the  aerodynamically  induced 
thermal  stresses  in  a  radome  constructed  of  Pyroceram  9606. 

The  analysis  was  to  be  conducted  on  ten  distinct  trajec¬ 
tories  with  thermal  stresses  calculated  at  three  separate  time 
steps  during  each  trajectory.  Temperature  dependent  material 
properties,  the  geometry  of  theoretical  radomes ,  and  the  thermal 
profiles  in  the  radomes  at  selected  time  steps  in  a  trajectory 
were  provided  by  the  AFML.  In  order  to  examine  many  geometric 
and  material  variations  as  well  as  the  thermal  conditions  re¬ 
sulting  from  various  aerodynamic  trajectories,  a  pre-processor 
was  written  which  used  the  geometric  and  temperature  output 
by  the  AFML  trajectory  analyses  as  input.  This  preprocessor 
generated  a  NASTRAN  input  deck  by  using  the  AFML  nodal  point 
temperatures  for  the  thermal  loads  while  using  the  same  nodal 
information  to  generate  the  NASTRAN  grid  points  and  element 
cards.  This  procedure  created  an  axisymmetric  finite  element 
model  using  the  axisymmetric  elements  in  NASTRAN,  primarily 
the  trapezoidal  CTRAPRG  and  CTRAPAX  elements,  although  some 
use  was  made  of  the  triangular  CTRIARG  and  CTRIAAX  elements. 

The  initial  analyses  consisted  of  radomes  subjected  to 
axisymmetric  aerodynamic  heating,  with  the  highest  temperature 
gradients  occurring  about  eight  seconds  after  launch.  Models 
using  both  the  CTRAPAX  and  the  CTRAPRG  elements  were  analyzed 
at  six,  eight  and  ten  seconds  after  launch.  These  results 
showed  high  stresses  at  the  radome  nosetip.  Evaluation  of  the 
thermal  fields  showed  higher  temperature  gradients  in  the  axial 
direction  than  in  the  thickness  direction  of  the  radome  due  to 
the  geometry  of  the  elements.  A  finer  AFML  model  was  therefore 
generated  which  would  cut  the  axial  gradient  to  about  300°F 
instead  of  the  600°F  which  existed  in  the  first  model.  This 
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reduced  the  stresses  some,  but  they  still  appeared  to  be  too 
high,  so  a  further  iteration  was  made  reducing  the  maximum 
temperature  between  any  two  elements  to  approximately  200°F. 

This  was  about  as  fine  as  the  AFML  thermal  model  could  go  since 
it  was  limited  to  around  380  nodes.  In  addition,  the  aspect 
ratios  in  the  area  of  interest  had  been  reduced  to  approximately 
one-to-one.  For  this  model,  the  CTRAPRG  solution  predicted 
maximum  compressive  thermal  stresses  of  about  21000  psi,  which 
agreed  well  with  a  SAAS  III  analysis  performed  by  the  AFML. 

The  CTRAPAX  elements,  however,  predicted  stresses  as  much  as 
100%  higher. 

The  reason  for  the  variation  in  the  stresses  predicted 
by  the  two  elements  was  unknown.  According  to  theory,  the 
stresses  should  have  been  identical  for  axisymmetric  loading. 
Various  ways  of  modeling  the  radome  with  triangular  elements 
were  also  investigated  to  determine  if  the  problem  was  a  func¬ 
tion  of  modeling  techniques.  This  did  not  appear  to  be  so  since 
all  combinations  of  AX  elements  produced  essentially  similar 
answers,  and  the  RG  elements  likewise  produced  a  set  of  similar 
stresses.  It  was  necessary  to  pursue  this  discrepancy  further 
since  the  AFML  requested  ASIAC  to  investigate  thermal  stresses 
due  to  non-axisymmetric  aerodynamic  heating  instead  of  contin¬ 
uing  the  symmetrically  loaded  thermal  cases.  Since  only  the 
CTRAPAX  and  CTRIAAX  elements  could  handle  asymmetric  loading, 
it  was  necessary  to  determine  the  reliability  of  these  elements. 
For  this  purpose,  a  simple  disc  model  with  known  theoretical 
solutions  was  examined  [1]".  A  quadratically  varying  radial 
temperature  and  a  linearly  varying  axial  temperature  were  ex¬ 
amined  to  see  if  the  discrepancy  was  due  to  the  two-dimensional 
nature  of  the  thermal  field  being  applied.  For  temperature 
independent  material  properties,  CTRAPRG  and  CTRAPAX  models 
agreed  exactly  with  the  theoretical  solution  as  long  as  all 


^Numbers  in  brackets  indicate  references  at  end  of  report. 
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elements  were  rectangularly  shaped.  When  the  shapes  were 
changed  to  a  non-rect angular  parallelogram,  the  CTRAPAX  pre¬ 
dicted  radial  stresses  60  to  80  percent  higher  than  the  theo¬ 
retical,  while  the  CTRAPRG  stresses  remained  exact.  Similar 
runs  with  temperature  dependent  material  properties  produced 
similar  results.  Even  worse,  the  CTRAPAX  model  predicted  axial 
stresses  as  high  as  27000  psi,  when  the  theoretical  solution 
and  the  CTRAPRG  results  predicted  no  axial  stress.  This  effect 
was  shown  to  be  essentially  aspect  ratio  independent ,  depending 
almost  entirely  on  the  degree  of  skewness  of  the  parallelograms. 

Since  it  was  very  difficult  to  model  the  radomes  without 
using  skewed  triangles  or  parallelograms,  a  different  approach 
was  needed.  The  cyclic  symmetry  option  in  NASTRAN  was  examined. 

A  ring  model  was  constructed  using  the  solid  CIHEXI  elements 
with  four  elements  in  the  axial  direction,  eight  in  the  radial 
direction,  and  one  in  the  circumferential  direction,  forming 
^  thirty  degree  wedge.  Theoretical  solutions  exist  for  a  ring 
subjected  to  a  temperature  distribution  of  the  form  Tq  r  cos 
n6,  where  k  and  n  are  integers  [2].  The  NASTRAN  results  were 
compared  to  the  theoretical  predictions  for  both  rectangular 
shaped  and  skewed  parallelepiped  elements.  These  produced 
almost  identical  stresses  so  this  technique  does  not  appear  to 
be  sensitive  to  the  geometry  used  for  the  elements.  The  stresses 
were  about  ten  percent  lower  than  theoretical ,  but  this  can  be 
attributed  to  the  large  arc  used  to  model  the  wedge.  Reducing 
the  wedge  to  ten  degrees  produced  much  better  results,  but  the 
computer  time  required  increased  proportionately. 

When  compared  to  the  axisymmetric  elements ,  the  cyclic 
symmetry  option  appears  to  offer  an  excellent  alternative  for 
the  solution  of  axisymmetric  structures  subjected  to  asymmetric 
heating  such  as  required  in  this  type  of  radome  analysis. 
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SECTION  2 
RADOME  ANALYSIS 

The  basic  radome  shape  used  in  the  analysis  is  shown  in 
Figure  1.  The  plot  is  approximately  actual  size,  except  that 
it  is  shortened  by  omitting  most  of  the  center  section.  The 
a^Pual  model  simulated  was  an  ogive  19.6  inches  long,  6.5  inches 
in  diameter  at  the  base  and  .281  inches  thick.  This  model  con¬ 
tained  9  nodes  through  the  thickness  and  40  in  the  axial  direc¬ 
tion,  with  a  total  of  339  grid  points  and  296  CTRAPRG  and 
CTRIARG  elements.  The  spacing  of  the  nodal  points  was  dictated 
by  the  AFML  program  which  calculated  the  temperature  distribution 
as  a  function  of  the  trajectory.  This  program  output  tempera¬ 
ture  at  the  centers  of  seven  elements  through  the  thickness, 
plus  the  temperature  at  the  inner  and  outer  surface  of  the 
shell.  These  points  were  then  used  as  NASTRAN  grid  points  with 
the  corresponding  nodal  point  temperatures. 

The  aerodynamic  heating  created  high  temperature  gradients 
in  both  the  axial  and  radial  direction.  Because  of  the  tempera¬ 
ture  dependent  material  properties,  it  was  desirable  to  keep 
these  gradients  as  small  as  possible.  Figure  2  represents  the 
best  distribution  of  thermal  gradients  that  was  obtainable  at 
the  nosetip.  This  distribution  produced  the  circumferential, 
radial,  axial,  and  shear  stresses  shown  in  Figures  3  through  6. 
These  figures  also  show  the  element  selection  that  produced  the 
most  consistent  set  of  stresses.  Comparison  with  Figure  2, 
which  was  the  original  element  orientation,  shows  that  the 
triangular  elements  have  been  moved  to  the  inner  surface  of  the 
model  instead  of  being  at  the  outer  surface.  It  can  be  observed 
•tha-t;  the  stresses  in  the  triangles  do  not  appear  consistent  with 
the  stresses  in  the  neighboring  trapezoidal  elements.  Several 
other  types  of  modeling  were  tried,  including  all  triangular 
elements;  however,  the  triangular  elements  always  appeared  to 
produce  stress  fields  with  sudden  changes  or  stress  reversals. 
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Figure  3.  Circumferential  Stresses,  psi,  in  Radome  Nosetip. 
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Figure  4.  Radial  Stresses,  psi,  in  Radome  Nosetip 
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Figure  6.  Shear  Stresses,  psi,  in  Radome  Nosetip 


CTRAPAX  and  CTRIAAX  elements  were  also  examined,  but  these 
elements  exhibited  stresses  as  high  as  30,000  to  40,000  psi  for 
the  same  temperature  distributions. 

The  AFML  accepted  the  results  as  shown  in  Figures  3,4,  5 
and  6  as  being  consistent  with  stresses  calculated  by  SAAS  III, 
and  thus,  decided  that  SAAS  III  could  be  used  to  complete  the 
analysis  of  the  axisymmetric  heating.  Based  on  the  stress 
levels  shown  by  this  analysis,  AFML  decided  that  the  nosetip 
would  be  removed  from  the  new  model  and  replaced  by  a  steel  plug 
to  relieve  the  thermal  stresses.  These  results  and  the  deci¬ 
sion  to  insert  a  metal  nosetip  parallels  closely  a  similar  ana¬ 
lysis  of  a  silicon  nitride  radome  subjected  to  aerodynamic 
heating  using  SAAS  II  [3]. 
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SECTION  3 
DISC  ANALYSIS 

In  order  to  determine  why  the  CTRAPAX  and  CTRAPRG  elements 
produced  differing  results  for  the  radome  model,  a  simple  disc 
model  restrained  at  the  outer  circumference  which  had  known 
theoretical  solutions  was  chosen  [1].  Both  a  quadratically 
varying  radial  temperature  and  a  linearly  varying  axial  tempera¬ 
ture  could  be  applied.  This  model  and  a  cross-section  of  the 
axisymmetric  model  is  shown  in  Figure  7(a).  For  this  model, 
both  types  of  elements  produced  exact  theoretical  answers  for 
aspect  ratios  varying  from  1.0  to  10.0  for  temperature  indepen¬ 
dent  material  properties  for  combinations  of  radially  and  axially 
varying  temperatures.  This  eliminated  aspect  ratios  and  two- 
dimensional  temperature  gradients  as  being  responsible  for  the 
differing  results  in  the  radome.  There  was  no  theoretical 
solution  for  temperature  dependent  material  properties;  however, 
both  elements  still  produced  identical  answers.  This  also 
eliminated  temperature  dependent  properties  as  being  the  cause, 
qg^ving  geometry  of  the  elements  as  the  only  likely  remaining 

candidate . 

The  geometry  of  the  elements  was  changed,  as  shown  in 
Figure  7(b),  so  that  all  elements  had  at  least  one  skewed  side 
instead  of  being  rectangular.  This  model  was  then  run  with 
only  the  radially  varying  temperature  distribution.  These  runs 
produced  results  that  definitely  proved  that  the  CTRAPAX  ele¬ 
ments  produced  incorrect  results  when  a  non-rectangular  cross- 
section  is  used.  For  the  temperature  independent  results,  the 
CTRAPRG  elements  produced  results  which  matched  the  theoretical 
solution  exactly,  while  the  CTRAPAX  elements  gave  radial  and 
hoop  stresses  sixty  to  one  hundred  percent  too  high.  Even 
worse,  these  elements  predicted  axial  stresses  almost  as  high 
as  the  axial  and  radial  stresses,  while  the  CTRAPRG  results 
agreed  with  the  theoretical  solution  of  zero  stress.  The 
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Figure  7(b).  Axisymmetric  Disc  Model  with  Skewed  Elements 
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temperature  dependent  material  runs  predicted  stresses  that 
followed  the  same  pattern  but,  of  course,  could  not  be  compared 
to  the  theoretical  solution.  These  results  are  shown  in  Figures 
8  through  13. 

This  analysis  showed  that  the  CTRAPRG  elements  appeared  to 
produce  reliable  results  while  confirming  that  the  CTRAPAX 
elements  could  not  be  trusted  in  a  model  requiring  the  use  of 
non-rectangular  element  shapes.  This  essentially  ruled  out  the 
use  of  the  CTRAPAX  elements  in  a  model  simulating  a  radome 
shape.  Since  only  the  CTRAPAX  and  CTRIAAX  elements  can  be  used 
for  axisymmetric  models  subjected  to  non-axisymmetric  loads, 
it  was  necessary  to  look  for  an  alternate  way  of  solving  the 
asymmetric  aerodynamic  heating  problem. 
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Ficure  13.  Axial  Stress  (Temperature  Dependent 

B  Materials)  for  DISC  with  Skewed  Elements 


SECTION  4 
RING  ANALYSIS 


The  cyclic  symmetry  option  in  NASTRAN  was  examined  to  de¬ 
termine  if  better  results  for  symmetric  structures  subjected  to 
an  asymmetric  thermal  load  could  be  obtained.  A  ring  subjected 
to  a  temperature  distribution  of  the  form  T  =  T^  r  cos  n0  was 
chosen  because  theoretical  solutions  could  be  obtained  [2]. 

The  model  of  the  ring  is  shown  in  Figure  14.  Both  thirty  degree 
and  ten  degree  wedge  shapes  were  examined,  requiring  twelve  and 
thirty-six  load  cases,  respectively,  when  running  cyclic  sym¬ 
metry.  The  model  cross  section  was  deliberately  made  as  similar 
to  the  previous  disc  analysis  as  possible,  including  the  use 
of  skewed  elements  shapes.  Table  1  shows  the  temperatures  used 
Input  and  the  resulting  theoretical  hoop  and  radial  stresses. 
The  axial  stress  should  be  identically  zero.  The  actual  stresses 
obtained  at  several  radii  are  shown  in  Tables  2  and  3  for  both 
ten  and  thirty  degree  wedges  with  skewed  and  non-skewed  geometry. 
It  can  be  seen  that  skewness  had  little  effect  except  for  the 
radial  stress  in  the  outermost  elements.  The  axial  stresses 
tended  to  be  less  than  ten  percent  of  the  lower  of  the  radial 
or  hoop  stress,  except  at  the  outer  fiber.  It  was  discovered 
that  the  axial  stresses  could  be  made  smaller  by  making  the 
ring  thinner;  thus,  approaching  a  state  of  plane  stress  more 
closely.  Selected  plots  of  hoop  stress  are  shown  as  Figures 
15,  16  and  17.  It  can  be  seen  in  these  figures  that  as  the 
wedge  becomes  narrower,  it  appears  to  approach  the  theoretical 
solution  as  a  limit. 

In  order  to  make  a  comparison  of  the  computer  costs  of 
cyclic  symmetry  against  the  use  of  the  axisymmetric  elements, 
a  ring  model  with  the  same  geometry  using  the  CTRAPAX  elements 
with  non-skewed  geometry  and  temperature  input  at  every  fifteen 
degrees,  as  shown  in  Table  1,  was  examined.  This  produced 
almost  exact  theoretical  answers;  however,  the  axisymmetric  model 


21 


T  = 

»-3 

o 

N5 

O 

o 

CO 

CD 

orr 

=  E  -  T0  Bk 

COS0 

a66 

=  E  cc  to  Dk 

COS0 

r 

Bk 

Dk 

72 

0.0 

.1202 

.3 

.03258 

.1177 

.4 

.04408 

.1163 

.  5 

. 04808 

.1024 

.  6 

.04696 

.07364 

.7 

. 04144 

.02929 

.8 

.03170 

-.03088 

.9 

.01788 

-.1070 

1.0 

0.0 

-.1990 

CIHEX  WEDGE  MODEL  USED  FOR  CYCLIC  SYMMETRY 
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Figure  14.  Axi symmetric  Ring  Model  Subjected  to 

Asymmetric  Heating 
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Table  2.  STRESSES  IN  30°  WEDGE  MODEL  AS  A 
FUNCTION  OF  ANGLE  AND  RADIUS 


Ring  With  Ring  With 

Un shewed  Elements  Shewed  Elements 


STRESSES 

AT  .25 

INCH 

RADIUS 

THETA 

(Degrees) 

RADIAL 

HOOP  AXIAL 

RADIAL 

HOOP 

AXIAL 

0 

41 

187 

-4 

42 

187 

-3 

30.0 

35 

162 

-4 

36 

162 

-2 

60.0 

20 

93 

-2 

21 

94 

-2 

90.0 

0 

0 

0 

0 

0 

0 

STRESSES 

AT  .55 

INCH 

RADIUS 

0 

87 

158 

8 

87 

158 

9 

30.0 

76 

137 

7 

76 

137 

8 

60.0 

43 

79 

4 

44 

79 

5 

90.0 

0 

0 

0 

0 

0 

0 

STRESSES 

AT  .95 

INCH 

RADIUS 

0 

26 

-221 

43 

26 

-221 

43 

30.0 

23 

-192 

37 

22 

-192 

37 

60.0 

13 

-111 

22 

13 

-111 

21 

90.0 

0 

0 

0 

0 

0 

0 
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Table  3.  STRESSES  IN  10°  WEDGE  MODEL  AS  A 
FUNCTION  OF  ANGLE  AND  RADIUS 


Ring  With  Ring  With 

Un shewed  Elements  Shewed  Elements 


STRESSES 

AT  . 

25  INCH 

RADIUS 

THETA 

RADIAL 

HOOP  AXIAL 

RADIAL 

HOOP 

AXIAL 

(Degrees) 

0.0 

46 

224 

-8 

50 

224 

-5 

30.0 

40 

194 

-7 

43 

194 

-5 

60-0 

23 

112 

-4 

25 

112 

-3 

90-0 

0 

0 

0 

0 

0 

0 

STRESSES 

AT  . 

55  INCH 

RADIUS 

0.0 

101 

178 

0 

102 

180 

2 

30.0 

87 

154 

0 

89 

156 

1 

60.0 

50 

89 

0 

51 

90 

1 

90.0 

0 

0 

0 

0 

0 

0 

STRESSES 

AT 

.95  INCH 

RADIUS 

0.0 

28 

-276 

33 

18 

-279 

33 

30.0 

25 

-239 

28 

16 

-242 

29 

60.  0 

14 

-138 

16 

9 

-140 

17 

90-0 

0 

0 

0 

0 

0 

0 

25 


THEORETICAL  STRESSES 
STRESSES  FROM  30°  WEDGE 
STRESSES  FROM  10°  WEDGE 


Stress  as  a  Function  of  Angle  for  a  Radius 
25  Inches  in  the  Asymmetrically  Heated  Ring 


THEORETICAL  STRESSES 
STRESSES  FROM  30°  WEDGE 
STRESSES  FROM  10°  WEDGE 
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Figure  16.  Hoop  Stress  as  a  Function  of  Angle  for  a  Radius 
of  .55  Inches  in  the  Asymmetrically  Heated  Ring 


required  that  the  CTRAPAX  elements  be  rectangular,  while  the 

symmetry  option  did  not  have  this  limitation.  The  fol¬ 
lowing  is  a  comparison  of  the  running  time  on  level  16.0  NASTRAN 
with  32  elements  in  each  model. 


CYCLIC  SYMMETRY  TECHNIQUE 

(30°  Wedge)  (10°  Wedge) 


AXI SYMMETRIC  TECHNIQUE 
(15°  Increments) 


CM (octal)  165,000  170,000 
CP (sec)  445  1,200 
10 ( sec )  333  869 


250,000 

2,681 

307 


These  results  indicate  that  the  cyclic  symmetry  option 
in  NASTRAN  is  better  suited  to  the  solution  of  a  general  axi- 
symmetric  problem  under  asymmetric  loading  than  the  axisymmetnc 
technique . 
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